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We propose and investigate an optical scheme for probabilistic implementation of an arbitrary 
single-mode quantum operation that can be expressed as a function of photon number operator. 
The scheme coherently combines multiple photon addition and subtraction and is feasible with 
current technology. As concrete examples, we demonstrate that the device can perform approximate 
noiseless linear amplification of light and can emulate Kerr nonlinearity. 
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I. INTRODUCTION 

Quantum properties of light have attracted the atten- 
tion of scientists since the early days of quantum me- 
chanics. Arguably one of the main goals in the field 
of quantum optics is to realize highly nonlinear inter- 
actions at the few-photon level which would enable e.g. 
to generate various highly nonclassical states of light and 
implement advanced schemes for quantum information 
processing. Unfortunately, nonlinear coupling between 
single photons mediated by common material media is 
extremely weak so other approaches have to be pursued. 
One of the most promising techniques appears to be that 
of the measurement induced nonlinearities. As shown by 
Knill, Laflamme, and Milburn using passive linear op- 
tics, ancilla single photons and photon counting measure- 
ments one can emulate nonlinear coupling between single 
photons and implement all-optical quantum CNOT gate. 
During recent years, the feasibility of this approach has 
been corroborated by numerous experiments 0. More- 
over, these ideas have been extended also to the so-called 
continuous variable regime, where quantum information 
is encoded into states of optical modes instead of single 
photons. It was shown that linear canonical transforma- 
tions of quadrature operators such as squeezing or quan- 
tum non-demolition coupling can be implemented with 
passive linear optics, ancilla squeezed vacuum states, ho- 
modyne detection and feedforward [H, 0, IE]- Going be- 
yond Gaussian operations, the basic techniques available 
arc the addition [6] and subtraction Q of a single pho- 
ton. These operations enable to generate highly non- 
classical superpositions of coherent states [1, UHo], and 
distill and concentrate continuous- variable entanglement 
[TTI . Il2| . If these operations are combined with coherent 
displacement, it is possible to prepare an arbitrary single- 
mode state of light from initial vacuum or squeezed state 

[SQ. 

As both single-photon subtraction and addition have 
been successfully demonstrated experimentally, it is in- 
teresting to investigate whether these elementary trans- 
formations can be combined to engineer some more com- 
plex useful quantum operations on the states of travel- 
ing light beams. In this paper we propose and analyze 
scheme for approximate probabilistic realization of an ar- 



bitrary operation that can be expressed as a function of 
photon number operator n. This class of transformations 
includes for instance the Kerr nonlinearity described by 
a unitary operation U = exp(— itfin 2 ), or a noiseless lin- 
ear amplifier [IE, HE, [TtJ Z = g", where g > 1 is the 
amplification gain. Our proposal is inspired by Ref. [IE] , 
where a scheme that coherently combines photon sub- 
traction and addition has been devised for the purpose 
of direct verification of the bosonic commutation rela- 
tions for creation and annihilation operators. We show 
that this scheme can be generalized to realize arbitrary 
operation f(n). An essential advantage of our approach 
is that we preserve the relative simplicity of the setup 
discussed in Ref. [IE]. The scheme involves only a sin- 
gle Mach-Zchnder interferometer and a single nonlinear 
crystal where parametric down-conversion occurs. The 
precision of the approximation of a given operation /(n) 
is controlled by the number of photons N counted by 
two photo-detection blocks that form a part of the setup. 
In contrast, other proposals for emulation of operations 
f(h) require either serial [l9j] or parallel [IE, [Ty] imple- 
mentation of many basic building blocks. Our approach 
is thus appealing from the experimental point of view, 
especially given that the core scheme has already been 
successfully realized experimentally by Zavatta, Bellini 
and coworkers [IE, Hi] . 

The rest of the present paper is organized as follows. 
In Sec. II we describe the proposed scheme and show how 
to determine its parameters for a given target operation 
f(n). In Sec. Ill we consider approximate realization 
of a noiseless amplification. In. Sec. IV we study the 
emulation of Kerr nonlinearity. Finally, conclusions are 
drawn in Sec. V. 



II. SETUP DESCRIPTION 

The scheme for probabilistic implementation of an ar- 
bitrary transformation f(n) that can be expressed as a 
function of photon number operator n is depicted in Fig. 
1. As already mentioned in the introduction, this scheme 
is very similar to the setup proposed in Ref. [Hj for ex- 
perimental testing of the bosonic commutation relation 
for annihilation and creation operators, [a, o^] = 1. Here 
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FIG. 1: (Color online) Linear optical implementation of oper- 
ators that are polynomials in photon number operator n. The 
scheme consists of unbalanced beam splitters BSi and BS2, 
polarizing beam splitter PBS, half- wave plate HWP, mirror 
M, nonlinear crystal NLC pumped by a strong laser pulse P, 
and two detection blocks Di and D2. The detection block Di 
counts the number of photons in mode C. Detection block 
D2 projects the state in mode B onto a specific entangled 
TV-photon polarization state. For more details, see text. 



squeezing operator, 
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Here a and c denote the annihilation operators of signal 
and idler modes, respectively, and A = tanhs where s is 
the squeezing constant. Assuming that the idler mode 
is initially in the vacuum state and that the detector Di 
registers N photons, the conditional transformation of 
the state of signal mode reads 
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In the limit of low parametric gain, A — > 0, we obtain 
An oc aJ N as desired. 

The photon subtraction is performed by splitting a tiny 
part of the signal beam off a highly unbalanced beam 
splitter followed by detection of the number of reflected 
photons @, H, [I] • If k photons are reflected from a beam 
splitter with amplitude transmittance t and reflectance 
r then the conditional transformation of the signal-mode 
state can be expressed as follows, 
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we consider a generalized version of that setup where N 
photons are added to and also N photons are subtracted 
from the input state. 

The photon addition is achieved by feeding the state 
into the input signal port of a nonlinear crystal NLC 
where pairs of correlated signal and idler photons are gen- 
erated in the process of noncollincar parametric down- 
conversion pumped by a strong coherent laser beam P 
The photodetector Di counts the number of emit- 
ted idler photons N which is equal to the number of 
photons added into the signal mode. Photon counting 
can be achieved e.g. by employing the time-multiplexed 
photon number resolving detector [22|, HH, [24[ where the 
input pulse is divided into a sequence of many pulses by 
highly unbalanced interferometers. The time-separated 
pulses then impinge onto avalanche photodiodcs (APD) 
and the total number of clicks of APDs gives the number 
of counted photons N. The main limitations of this detec- 
tion scheme are the finite number of time bins which does 
not allow for complete photon number resolution, and 
the non-unit detection efficiency of the APDs. However, 
these limitations are not very restrictive for the present 
application where the down-conversion should be oper- 
ated in a regime of low gain. The probability that N + 1 
photons are present in the output idler beam is then much 
lower than the probability that N photons are present. If 
N photons are counted by Di then, with high probabil- 
ity, N photons were indeed emitted into the idler mode 
in NLC. In this regime, the main effect of the less-than- 
unit detection efficiency of APDs is a reduction of success 
probability of the scheme. 

Mathematically, the state transformation after propa- 
gation through the crystal is described by the two-mode 



A crucial feature of the scheme in Fig. 1 is that the 
photon subtraction may occur either before or after the 
photon addition. The photons reflected off beam splitters 
BSi and BS2 are recombined on a polarizing beam split- 
ter PBS such that their polarization states are orthogonal 
but all other degrees of freedom are made indistinguish- 
able to guarantee maximum visibility of multiphoton in- 
terference. The detector D2 represents a detection block 
that is capable of making projection onto an TV-photon 
polarization state 



jv 



)b =J2b k \k) BtH \N -k) B ,v, 



(4) 



k=0 



where H and V denote the horizontally and vertically 
linearly polarized modes, respectively. Such a projec- 
tion can be performed e.g. by reversing a linear optical 
scheme for preparation of the A^-photon two-mode states 
put forward in Ref. [2||. The input light beam is di- 
vided into A" spatial modes and a polarization analysis 
is performed on each mode by an elementary detection 
block consisting of quarter- and half- wave plates, polariz- 
ing beam splitter and a pair of single photon detectors. If 
a single photon is detected by each detection block, and 
a photon in the j-th detection block is projected onto a 
polarization state 



cos 9j \H) B + sin 6j e 1 ^ \V) B , 



(5) 



then the whole detector projects onto A^-photon polar- 
ization state [25|, [26[ 
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Here b\j, b v denote creation operators of the horizontally 
and vertically polarized modes, respectively. Since every 
A-photon polarization state can be expressed in the fac- 
torized form ([5]), this detection scheme can project onto 
any state (|4]). 

If N photons in total are subtracted by BSi and BS2 
and if addition of N photons is heralded by the detector 
Di, then the operation on signal mode A and the N- 
photon polarization state of auxiliary spatial mode B are 
given by, 
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With a slight abuse of notation, we use in Eq. (O the 
tensor product symbol ® to separate operator acting on 
the Hilbert space of the signal mode and the A-photon 
state prepared in the auxiliary mode B. After projecting 
the mode B onto the N-photon state (j4]) we finally obtain 
the operation on signal mode, 
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where C N = VI - A 2 (rX) N /N\ and h = tVl - A 2 . 
With the help of the expression 
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we can rewrite Wn as follows, 

W N = C N h A P N (n), (10) 
where P/v(n) is a polynomial of A-th order in n, 
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We now prove that by properly tailoring the state ([!]), 
arbitrary polynomial P/v(n) can be obtained. We start 
from Eq. ([5]) where we normally order the creation and 
annihilation operators using the formula, 



3=0 

In this way we obtain 
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where dj = Y, k =i Mjkbk, and 
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By setting Mjk = 0, j > fc, we can express the relation 



between dj and b^ in a matrix form, 



d = Mb, 



(15) 



where b = (bo, bi, ... , 6at) t , d = (do, d%, . . . , c?at) t . The 



system of equations (|T5j) for bk can be always solved. 
Since Mjj > 0, Vj and Mjk = 0, j > k, we have 
dct M > 0, M^ 1 exists and we can write b = M d. 
Any polynomial in h can be recast into the form (|13p by 
normally ordering the annihilation and creation opera- 
tors and the coefficients d k thus unambiguously specify 
the polynomial. The amplitudes bk yielding the required 
dj can be calculated as described above which completes 
the engineering of the operator polynomial P/y(n). 

To summarize our findings so far, we have shown that 
with the scheme shown in Fig. 1 we can implement an 
arbitrary operator that can be expressed as a product of 
a polynomial in n and an attenuation factor h n (recall 
that h < 1). The scheme can approximate any given op- 
erator f(h) to arbitrarily high degree N determined by 
the number of photon subtractions and additions. More 
specifically, the first N + 1 terms in Taylor series expan- 
sion in h of f(n) and Wm can be made equal up to a con- 
stant prefactor. However, the success probability of the 
scheme will decrease approximately exponentially with 
increasing N . Practical implementations of the proposed 
scheme in a near future would be thus most probably 
limited to N = 1 and N = 2. 

In this low-approximation regime the influence of the 
attenuation factor h n could be particularly significant. 
This effective attenuation can be suppressed by employ- 
ing highly unbalanced beam splitters, t — ► 1, and work- 
ing in the regime of very low parametric gain, A — ► 0. 
Then we get Wn — Cn Pn ("-)■ However, in this limit 
the success probability of the scheme becomes very low 
and vanishes when t = 1 or A = 0. In an experiment, 
optimal balance between success rate and performance 
of the scheme can be sought by tuning A and t. 



III. NOISELESS AMPLIFIER 

As a first application of the scheme let us consider 
probabilistic implementation of a noiseless linear ampli- 
fier El El, 



Z(g) = /, 



(16) 



where g > 1 denotes an amplitude gain of the ampli- 
fier. The amplifier (|16p acts as a non-unitary quantum 
filter that enhances amplitudes of the Fock states \n) by 
a factor g n . Clearly, this transformation cannot be im- 
plemented exactly with a finite probability because the 
operator Z(g) is unbounded for any g > 1. We can never- 
theless approximately realize the transformation to A-th 
order in n by the scheme described in a previous section. 
In essence in this approach the amplification works well 
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FIG. 2: Dependence of the effective gain g e g of the prob- 
abilistic amplifier on the amplitude of coherent state \a\ is 
plotted for five different nominal values of the gain, g = 
1.1, 1.2, 1.3, 1.4, 1.5 and for two levels of approximation N = 1 
(a) and N = 2 (b). 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 




0.99 
0.988 
0.986 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

M 

FIG. 3: The fidelity of the amplified coherent states is plot- 
ted as a function of the coherent state amplitude a for five 
different gains g = 1.1, 1.2, 1.3, 1.4, 1.5 and for N = 1 (a) and 
TV = 2(b). 



It is natural to require that g e s = g in the limit \a\ — > 0. 
This can be achieved if we set d to be equal to a root of 
a polynomial equation 



for low Fock states \n), n < N, but it fails for large Fock 
states n^> N. 

For the sake of simplicity, in what follows we shall as- 
sume that h = 1, i.e. the limit of a low reflectance of 
beam splitters BSi, BS2 and low gain of the nonlinear 
crystal. The influence of h < 1 on the performance of 
the amplifier is that it reduces the effective gain, g — > hg. 
This can be compensated by substitution g — > g/h when 
determining the parameters of the scheme. 

We can construct a polynomial approximation to the 
operator Z{g) by expanding it in Taylor series in n and 
keeping the first N + 1 terms. In this way we obtain a 
polynomial 
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where d = hig. An important feature of the amplifier 
(fPo)) is that it preserves the structure of coherent states 
and only amplifies their amplitude, Z\a) oc \ga). It turns 
out that the gain of the approximate amplifier (fT"7]) de- 
pends on I a I and we can define an effective gain of the 
amplifier, 
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As a concrete example we now study in more detail the 
cases N = 1 and N = 2 that are experimentally feasible 
with current technology. For N — 1 we find that d = g— 1 
and the approximate polynomial reads 



Z x = (g-l)n+l. 



(20) 



For N = 2we must solve a quadratic equation for d which 
yields d — y/2g — 1 — 1 and 



(g - v/2.9 - l)n 2 + (^2g - 1 - l)n + 1. (21) 
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On inserting the operators ((20)) and (f2Tj) into Eq. (fl~8| 
we can derive analytical formulas for the effective gain 
For N = 1 we obtain 



5off = 1 
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The expression for N = 2 is rather long and unwieldy 
and is not reproduced here. Instead, we plot in Fig. 2 
the dependence of g e s on \a\ for N = 1, N = 2 and for 
five different nominal values of the gain. We can see that 
the effective gain decreases with increasing \a\. Also, 
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FIG. 4: The figure shows dependence of fidelity (a) and en- 
tropy of entanglement (b) of the amplified two-mode squeezed 
vacuum state on the amplification gain. The parameters read 
£ = §, ft -> 1, N = 1 (solid line), and N = 2 (dashed line). 

the performance of the scheme improves when higher- 
order approximation is employed, and we find that for 
N = 2 the effective gain g c g is closer to g than for N = 1 . 
Besides the effective gain, we also calculate the fidelity of 
the amplified weak coherent state with the ideal target 
state, 

p = \(ga\Z N \a)\ 2 
(a\Z f N Z N \a) 

We plot the fidelity as a function of coherent state am- 
plitude in Fig. 3. For the considered range of parameters 
g < 1.5 and \a\ > 1 the fidelity exceeds 0.95 for N = 1 
and surpasses 0.986 for N = 2. The scheme is thus very 
suitable for probabilistic error-free quantum cloning of 
weak coherent states via noiseless amplification. 

Another important application of a noiseless amplifier 
is the concentration of continuous variable entanglement 
[l5| . Consider a two- mode squeezed vacuum state 

oo 

\V)ab = ^Cn\n) A \n) B , (24) 

n=0 

where c n = \J\ — £ 2 £" and £ is the two-mode squeez- 
ing constant. Local noiseless amplification of one mode 
increases the two-mode squeezing, £ — > g£, and thus en- 
hances the entanglement of the state. We have investi- 
gated the effect of the approximate amplifications ([20]) 
and (f2"Tj) on the entangled state (f2"4"|) by means of nu- 
merical calculations. In the simulations we assume the 
limit h — ► 1. Generally, the attenuation factor h n can 
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FIG. 5: (Color online) The figure shows change of the Fock- 
state probability amplitudes c„ by noiseless amplification. For 
each n, the four bars from the left to the right correspond 
to the following states: initial two-mode squeezed state with 
squeezing constant £ (blue), state amplified by operation Z\ 
(red), state amplified by operation z?2 (green), target two- 
mode squeezed state with squeezing constant g£ (black). The 
parameters are £ = 3 and g = 1.75. 

be accounted for by replacing £ with rescaled squeezing 
constant £' = h£. 

The results are presented in Fig. 4 that shows the de- 
pendence of fidelity and entanglement of the amplified 
state on gain g for £ = •! which corresponds to 3 dB 
of two-mode squeezing. We can see that the amplifica- 
tion enhances the state entanglement as expected, while 
preserving well the structure of the state as witnessed by 
high fidelity. This is further illustrated in Fig. 5 where we 
plot the Fock-state probability amplitudes c n of a state 
locally amplified by operations Z\ and Z2 with g — 1.75, 
as well as the amplitudes of the initial two-mode squeezed 
vacuum state and the ideal amplified two-mode squeezed 
vacuum state with squeezing constant g£. 

IV. KERR NONLINEARITY 

The Kerr effect is a consequence of a nonlinear response 
of a medium whose index of refraction depends on the 
intensity of the light field. Quantum mechanically, this 
can be described by a Hamiltonian proportional to the 
square of photon number operator. The resulting unitary 
transformation reads, 

L> = e^™ 2 , (25) 

where is a dimensionlcss nonlinear phase shift. Kerr 
nonlinearity can be used to generate highly non-classical 
states of light beams and, together with Gaussian oper- 
ations, it is sufficient for universal quantum computing 
[27j |. With the scheme described in Sec. II we can imple- 
ment an approximate truncated version of the transfor- 
mation (f2"5"|) . 

fe=0 
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FIG. 6: (Color online) Plotted are the nonlinear phase shift 
$ n (blue triangles) and amplitude modulation factors A n 
(green squares) induced by the transformation (|28[l emulating 
Kerr nonlinearity with c/> = Also shown for comparison 
is the phase shift <jm 2 corresponding to the true Kerr nonlin- 
earity (|25)l (red circles). 



We can generalize this expression to compensate for the 
attenuation factor h n up to 2iV-th order by taking trun- 
cated Taylor series of operator e l ^ n h~ n , 
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where H^{y) denotes Hermite polynomial. In what fol- 
lows we shall assume the limit h — > 1 for the sake of sim- 
plicity. Emulation of Kerr nonlinearity is more resource- 
demanding than noiseless amplification, because already 
the first non-trivial approximation requires subtraction 
and addition of two photons, 



U 2 = 1 + idmr 



(28) 



The unitary matrix U is diagonal in Fock state basis, with 

'At 2 

matrix elements U nn = e tlpn . The matrix representing 
the approximate operation (j28|) is also diagonal in Fock 
state basis but we find that {U2) nn = \A + c/> 2 n 4 e 1 *™, 
where $„ = arctan(c/>n 2 ). The nonlinear phase shift $„ 
does not increase quadratically with n as it should but 
becomes saturated, $ n ^oo = f • Another aspect of the 
approximation (f2"8")) is that the Fock state amplitudes are 
modulated by the factors A n = yl + <j) 2 n 4 so the trans- 
formation acts also as an amplifier. These features are 
illustrated in Fig. 6. 

We quantify the performance of the approximate op- 
eration l[28p by the quantum process fidelity. Since the 
Hilbert space of a field mode is infinite dimensional, we 
shall consider fidelity Fn of operation restricted to a fi- 
nite dimensional subspace spanned by the first N+l Fock 
states |0), . . . , \ N). We define maximally entangled state 
on this subspace \^n) — ,^ J2 n =o l n )l n ) anc ^ we have 
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FIG. 7: Dependence of the quantum process fidelity Fn on 
the nonlinear phase shift cj> is plotted for N = 2 (solid line), 
N = 3 (dashed line) and N = 4 (dotted line). 



where 1 denotes the identity operation. Figure 7 shows 
the dependence of the fidelities F 2 , F$ and F4 on the 
phase shift tp. We can see that with growing (j> it is in- 
creasingly more and more difficult to emulate U by U 2 
and the fidelities Fn decrease rapidly. 

This problem of fidelity decrease can be remedied if 
we need to emulate the Kerr nonlinearity only on a finite 
dimensional subspace spanned by the first N + l Fock 
states. More generally, any operation f(n) can be per- 
fectly probabilistically implemented on such subspace if 
we construct the polynomial approximation to this oper- 
ation as follows, 
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Note that this expression includes compensation for the 
attenuation factor h. In particular, on inserting f(k) = 
e i0fc 2 and N _ 2 i nt0 Eq. fl3Tj]) we obtain the following 
polynomial that perfectly emulates Kerr nonlinearity on 
a three-dimensional subspace spanned by Fock states |0), 
|1), and |2), 
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(31) 



DISCUSSION AND CONCLUSIONS 



An important practical characteristics of the proposed 
device that we have not addressed in detail yet is its 
success probability. A simple order-of-magnitude esti- 
mate reads p S ucc ~ (rA?7) 2Ar , where 77 denotes the de- 
tection efficiency of single-photon detectors. The quan- 
tum efficiencies of APDs exceed 50%. However, in the 
experiment, the optical beams typically need to be spa- 
tially and spectrally filtered before detection which re- 
duces the overall detection efficiency to 10% or even less. 
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Another important characteristics is the attenuation fac- 
tor h = t 2 \/l — A 2 . We want to maximize p S ucc while 
keeping h as close to 1 as possible. To give a concrete 
example, we can set i 2 = 0.95, A 2 = 0.05, and r\ = 0.1. 
For N = 1 we obtain p succ » 2.5 x 10 -5 and h = 0.926. 
Pumping the scheme with a pulsed Tksapphire laser with 
a repetition rate of 80 MHz would then yield approx- 
imately 2000 events per second. Second-order approx- 
imation, N — 2, is much more challenging. For the 
above parameters we obtain success rate 0.05 Hz, i.e. 
about 3 events per minute which may be still accept- 
able for proof-of-principle experiments. Moreover, this 
rate can be significantly increased either by improving 
the detection efficiency or by using beam splitters with 
higher transmittance r. For instance, with t 2 = 0.9 and 
r\ = 0.2 we obtain success rate 3 Hz and h — 0.877. We 
stress that the calculated success rates are only order- 
of-magnitude estimates and would depend on the input 
state of the scheme as well as on details of the imple- 
mented transformation. Nevertheless, these calculations 
clearly confirm the experimental feasibility of the scheme 
with single photon addition and subtraction and show 
that scheme with two photon additions and subtractions 
may be within the reach of present technology. 
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